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Abstract 

The fundamental problem is analized, the relation between Bose-Einstein conden- 
sation and spontaneous gauge symmetry breaking. This relation is largerly misunder- 
stood in physics community. Numerous articles and books contain the statement that, 
though gauge symmetry breaking helps for describing Bose-Einstein condensation, but 
the latter, in principle, does not require any symmetry breaking. This, however, is 
not correct. The analysis is based on the known mathematical theorems. But in or- 
der not to overcomplicate the presentation and to make it accessible to all readers, 
technical details are often omitted here. The emphasis is made on the following basic 
general facts: Spontaneous breaking of gauge symmetry is the necessary and sufficient 
condition for Bose-Einstein condensation. Condensate fluctuations, in thermodynamic 
limit, are negligible. Their catastrophic behavior can arise only as a result of incor- 
rect calculations, when a Bose-condensed system is described without gauge symmetry 
breaking. It is crucially important to employ the representative statistical ensembles 
equipped with all conditions that are necessary for a unique and mathematically cor- 
rect description of the given statistical system. Only then one is able to develop a 
self-consistent theory, free of paradoxes. 
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1 Introduction 



In recent years, the phenomenon of Bose-Einstein condensation (BEC) has received much 
attention, both experimentally and theoretically (see, e.g., the book [1] and review arti- 
cles [2-6]). At the present time, Bose-Einstein condensates of trapped atoms have been 
realized in about 70 laboratories of 15 countries, USA, Germany, France, New Zealand, Eng- 
land, Japan, Italy, Netherlands, Israel, Australia, China, Austria, Switzerland, Canada, and 
Brazil. Magnetic, magneto-optical, and all-optical traps have been employed. BEC has been 
achieved for 11 atomic species (^H, ^He, ^Li, ^SNa, ^^K, ^^K, ^^Ci, ^^Rb, ^^Rb, ^^scs, and 
^^'''Yb), four types of molecules formed by Bose atoms (^^Na2, *^Rb2, ^^Rb2, and ^^^Cs2), 
and for two types of molecules formed by Fermi atoms (^Li2 and ^°K2). More details can be 
found in Refs. [2,5,7-10]. 

One of course, always remembers that BEC is assumed to exist in superfluid *^He, which 
can be revealed, e.g., in scattering experiments with neutrons [11-14] and x-rays [15,16]. 
Also, BEC of boson- type quark clusters might arise in interiors of neutron stars and in 
heavy-ion collisions [17-19]. There exist several other examples of BEC, which will be cited 
in the last section of the paper. Thus, Bose-condensed systems form quite a numerous family. 

In the theory of Bose-condensed systems, there have remained some principal problems 
that have not been properly understood in a large part of physics community. One of such 
fundamental problems is the relation between BEC and the U(l) gauge symmetry breaking. 
To understand the issue, it is necessary to give answers to the following group of mutually 
interrelated questions: Is gauge symmetry breaking only a sufficient condition for BEC or 
it is also a necessary one? In what sense one should interpret the statement that BEC and 
gauge symmetry breaking are equivalent? What is the difference between the two methods 
of gauge symmetry breaking, with the help of infinitesimal sources and by means of the 
Bogolubov shift? What is the origin of the catastrophic behavior of condensate fluctuations 
in Bose-condensed systems? And what arc the requirements for making a theory, with broken 
gauge symmetry, self- consistent? The aim of the present paper is to answer these questions. 

2 Expansions over natural orbitals 

We shall consider a statistical system of N bosons in volume V. For simplicity, we keep 
in mind the particles with no internal degrees of freedom. After answers to the principal 
questions are understood, it is easy to generalize the consideration to particles with internal 
degrees of freedom, when fragmented condensates [20] might arise. 

The basic variable, characterizing particles of a statistical system, is the field operator 
'0(r, t), where r is spatial vector and t, time. Bose particles correspond to the Bose commu- 
tation relations 



with other commutation relations for and [■0^,'0^] being zeros. 

2.1 Expansion of density matrix 

A general way for characterizing BEC is the Penrose- Onsager scheme [21], based on the 
consideration of the eigenvalues of the first-order density matrix, which is defined as the 




(1) 
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statistical average 

p(r,r',i) =<ij\r',t)ij{r,t) > . (2) 
The eigenvalues Uk are given by the eigenproblem 

j P(r, r', t) (pk{r', t) dr' = Ukit) (pk{r, t) , (3) 

whose eigenfunctions (fk are termed natural orbitals [22]. The density matrix (2) acquires 
the spectral resolution 

p{r,r',t) = Y.Mt) M^,t) ipl{r',t) (4) 

k 

over the natural orbitals. 

For an equilibrium system, the density matrix (2) does not depend on time. Respectively, 
Uk and (pk do not depend on time as well, so that the spectral resolution (4) becomes 

r') = ^k^k{r)^lir') . (5) 

k 

In Eqs. (3), (4), and (5), the index k is an appropriate multi-index labelling the eigenvalues 
and eigenfunctions of eigenproblem (3). 

In what follows, we shall mainly deal with the equilibrium case (5), keeping in mind that 
it can straightforwardly be extended to the general nonequilibrium situation of Eq. (4). 

Among the eigenvalues n^, let is separate the largest one 

No = supnk, (6) 

k 

related to a state (/?o(r) = V^fco('') labelled by a multi-index /oq. One says that there occurs 
BEG, when the largest eigenvalue Nq is proportional to the total number of particles N, so 
that Nq (X N for all N, including the thermodynamic limit, when 

N 

A^— >oo, y— >oo, — ^ const , (7) 

where a positive constant is implied. Thence, BEG, by definition, is characterized by the 
condition 

iVn 

Urn ^>0, (8) 

where the limit means the thermodynamic limit (7). Under this condition, the spectral 
resolution (5) can be rewritten as the sum 

p(r, r') = No<po{r)<p*o{r') + ^ nk<Pk{r)<pl{r') . (9) 

k^O 

The natural orbitals are assumed to be normalized to one, \\ipk\ \ — 1- Therefore, Nq is the 
number of condensed particles. 

Often one connects the appearance of BEG with the off-diagonal long-range order, defin- 
ing the condensate density as the limit 

Po = ^ = lim p(r, 0) . 
3 



This, however, has sense only for an equihbrium uniform system, when cpo = l/y/V and the 
second term in Eq. (9) tends to zero with r = |r| ^ oo. For example, for a nonuniform 
confined system, one has (/?fc(r) ^ as r — >■ oo, which would mean that for this nonuniform 
system BEC would be impossible [2]. 

The Penrose-Onsager criterion of BEC, in the form of condition (8), where A^o is the 
largest eigenvalue (6), is, clearly, more general. Another general form of the BEC criterion 
can be given by introducing the order indices [22,23]. Let us consider pi = [p(r, r')] as a 
matrix with respect to the variables r and r'. Then the order index of the matrix is defined 
[22,23] as 

N _ log||pi|| 

where ||pi|| is the norm of pi. For large ^ 1, this takes the form 

logiVo 
= 1^ ■ 

The order indices (10) or (11) are defined for an arbitrary system and describe different types 
of ordering: 

uj{pi) < (no order) , 

uj{pi) = (short — range order) , 

< uj{pi) < 1 (mid — range order) , 

u!{pi) — 1 (long — range order) . 

The long-range order happens in thermodynamic limit, when — > oo and A^o N. Then, 
from Eq. (11) it follows that 



hm a;(pi) = 1 . (12) 

AT— >oo 

The criteria (8) and (12) are, of course, absolutely general and can be applied to any 
system. However, being defined through the density matrix (2) , they are not useful without 
knowing the latter. But the calculation of the density matrix for an interacting system is 
not a trivial task. A practical calculational procedure is to be developed. 



2.2 Expansion of field operators 

Let us specify how the spectral resolution (9) for the density matrix can be obtained. The 
field operator ■0(r) can be expanded over the basis of natural orbitals, 

ip{r) = Y,ak^k{r) . (13) 

k 

The operators ak in the /c-representation obey the commutation relations 

[ofe, aj] = Skp , [ofe. Op] = . (14) 

Generally, we do not know which of the orbitals v2fc(r) corresponds to the condensate. 
The answer to this question is easy only for the uniform system, when the multi-index k 
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translates into the momentum k, so that the condensate is related to k = 0. But in the 
general case, we cannot decide which of the orbitals is that of the condensate. To decide 
this, we need to know the eigenvalues of the density matrix (2), with the largest eigenvalue 
pointing at the condensate. So that, actually, we are yet in the vicious circle. Expansion 
(13), however, is useful for proving some theorems of the condensate existence. 

Assume that we are aware which of the natural orbitals </?fe(r) is related to the condesate, 
and let us mark it as <^o(r). Then we may separate in expansion (13) the condensate field 
operator 

iJo{r) = ao9?o(r) (15) 
from the remaining part associated with the field operator of uncondensed particles, 

il^iir) = J2ak'Pk(r) . (16) 

k^O 

Hence, expansion (13) can be represented as the sum 

V^(r)=^o(r)+^i(r). (17) 
The natural orbitals are supposed to be orthonormal, 

/ fli^Wpi^) dr ^ 5kp . 

From here, the field operators of condensed and uncondensed particles, by construction, are 
orthogonal, 

' 4{r)ilj,{r) dr = Q . (18) 



I' 



Correctly speaking, the operators ipo{^) and V'i(r) are not separate independent operators, 
but they are just two parts of one field operator ip{r). All of these operators act on the Fock 
space J-{ip) generated by the operator ip"^ . All operators from the algebra of observables, 
being functionals of -0, are defined on the Fock space See mathematical details in Ref. 

[24]. 

Neither of the operators ipQ{r) or ^i(r) represents bosons, since their commutation rela- 
tions are not of Bose statistics. Really, the former operators obey the commutation relation 

[V'o(r), ^l:l{r')]^^o{v)vl{r'), (19) 

while the latter, the relation 

[V'i(r), V'I(r')] = E^^(rK(r'), (20) 

other commutation relations, such as [■0o,'0o], [V'l^V'i]) ^iid [■0o,'0i] being zeros. Involving 
the 5-function expansion 

5(r-r') = E^^(r)^:(r'), 

k 

that is, the completeness of the natural orbital basis {</7fe(r)}, the commutation relation (20) 
can be represented as 

[^i(r), ^l(r')] = 5{r - r') - [^o(r), ^5(r')] = 5{r - r') - ^o(r)^S(r') • (21) 
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Prom the latter equation it follows 

jlMr), ^I(r')]V^o(r') dr' = , 

in agreement with the orthogonality condition (18). 

The commutation relations (19), (20), and (21) simplify in thermodynamic limit, provided 
that the condensate natural orbital satisfies the inequality 

Mr)r<^ (.>0). (22) 

For illustration, we may recall that for a uniform system 

hence inequality (22) is valid, with i/ — 1. Another example is the gas of atoms trapped in 
a potential U{r) oc r". Then [25], inequality (22) is again vahd, with u — 1 + 2/n. Under 
condition (22), one has 

T I / N */ /M ^ 1- const 

hm ipo{r)(pQ{r) < lim — — = . 

AT— ►oo AT— >oo iV 

Therefore, in thermodynamic limit, one can formally write down Eq. (19) as 

lim[V^o(r), 4ir')] = 0, (23) 

N—*oo 

and the commutation relation (20) becomes of the Bose kind, 

lim[V^i(r), 4{r')] = d{r-r'). (24) 

N—*oo 

More rigorously, one should consider Eqs. (19) and (21) as equations for distributions, which 
have sense being integrated with integrable functions from the class of functions such that 



lim 

AT— »oo 



V 



/(r,r') drdr' 



< oo 



On this class of functions, taking into account condition (22), one gets 

hm /■ [V'o(r), 4ir^')]f{r,r')drdr'^0, 

AT— »oo Jv 

lim / [^i(r), xj;l(r')]f(r,r')drdr' ^ lim / S(r - r')f(r,r') drdr' , 

AT— »oo Jv AT— >oo Jv 

which explains Eqs. (23) and (24). 

Thus, in thermodynamic limit, the condensate variable (15) can be treated as commuting 
and the operators of uncondensed particles (16) as usual Bose field operators. But for a finite 
system, one has to deal with the exact commutation relations (19), (20), and (21). 

One says that BEG occurs in thermodynamic limit, when the number of condensed 
particles 

No^ J < '0o(r)V'o(r) > dr ^ < ajoo > (25) 



satisfies the condensation criteria (8) or (12). 

Substituting expansion (13), with notations (15), (16), and (17), into the density matrix 
(2), and assuming the quantum-number conservation condition 

< alup > = Skp< alttk > , 

we obtain the spectral resolution (9), in which A^o is given by Eq. (25), and 

rik ^ < alak > . 

Nowhere throughout this section, gauge symmetry breaking has been either invoked 
or even mentioned. The spectral resolution (9) formally does not require the breaking of 
symmetry [26]. The BEC criteria (8) and (12) also arc formulated without involving the 
notion of symmetry breaking. Because of this, one often states that BEC does not necessarily 
imply gauge symmetry breaking and that the overall phenomenon of BEC can perfectly be 
described without breaking the gauge symmetry. The latter conclusion, however, is not 
correct, since it does not follow from the fact that the spectral resolution (9) and the BEC 
criteria (8) and (12) can be formulated not mentioning symmetry breaking. Strictly specking, 
Eqs. (9), (8), and (12) tell us nothing about whether gauge symmetry must be broken or 
not. To understand this, it is necessary to know the correct way of calculating the density 
matrix and all other averages. 



3 Methods of symmetry breaking 

For the practical purpose of describing BEC, one usually breaks gauge symmetry. There 
are two ways of doing this, by the method of infinitesimal sources and by the Bogolubov 
operator shift. 



3.1 Method of infinitesimal sources 



Let the system Hamiltonian H[iIj], being a functional of the field operator ip, be invariant 
under the gauge transformation ip e*"-;/', where a is a real number, so that 



H 



(26) 



To break the gauge symmetry, we can add to the Hamiltonian H[ip] a term explicitly breaking 
the symmetry, for instance, defining 



H,[i;] = H[^]+e^ / ^^(r)+V^o(r 



dr 



(27) 



Statistical averages for an equilibrium system are then averaged with the statistical operator 

(28) 



. _ exp{-(3H,m 

Pe = 



in which /3 = 1/T is inverse tepmerature. Thus, the statistical average of an operator A is 
given by 

<A>,= TipeA, (29) 
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where the trace is over the Fock space ^(V")- 

Since Hamiltonian (27) is not invariant under the gauge transformation, the system 
symmetry is broken. The manifestation of the broken gauge symmetry is 

< V(r) >e^0. (30) 

But it is important to find out what happens if the symmetry breaking term is removed. 
If it is removed before thermodynamic hmit, then evidently 

hm < V'(r) >£ = . (31) 

However, the average of i/j may remain nonzero, if the symmetry breaking term is removed 
after thermodynamic hmit. One says that gauge symmetry is spontaneously broken, when 

hm hm < '^(r) >e ^ ■ (32) 

£— >0 W— »oo 

For a nonequihbrium system, the statistical operator (28) can be treated as the initial form 
of the time-dependent statistical operator Pe(t), so that ^^(O) = p^. 

The averages of type (32), defined through the double noncommuting limiting procedure 
limg^o liKiAf^oo < ••• >! arc the Bogohibov quasiaverages [27,28]. It is also possible to 
define the thermodynamic quasiaverages [29-31] involving the sole limiting procedure, that 
of thermodynamic limit lim^v^oo < ■ ■ ■ >■ The latter procedure can also result in the 
spontaneous symmetry breaking, provided the infinitesimal source in thermodynamic limit 
tends to zero in the appropriate way [29-31]. 



3.2 Bogolubov operator shift 

Another method of gauge symmetry breaking is based on the observation that in thermo- 
dynamic limit we come to the commutation relations (23) and (24). The first tells us that 
the operator ipo{r) is asymptotically equivalent to a nonoperator function, which can be 
expressed as the replacement 

V'o(r) ^ 77(r) , ao ^ , (33) 

where ri{r) is termed the condensate wave function. Respectively, the field operator (17) is 
to be replaced as 

^{r) ^ ^(r) = 77(r)+V'i(r), (34) 

which is called the Bogolubov operator shift. The commutation relation (24) signifies that 
the operator ^i(r) asymptotically acquires the usual Bose commutation relations. This 
agrees with the fact that if ip{r) is a Bose operator and rj^r) is a nonoperator function, then 
ipi{r) in Eq. (34) has to be a Bose operator with the commutation relation 

[Mr), ^^^(r')]=5(r-r'). (35) 
The orthogonality condition (18) transforms into 

J 77*(r)V'i(r) dr^O. (36) 
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Thus, after the Bogolubov shift (34), instead of one field operator ip{r), as in Eq. (17), there 
appear two independent variables, the condensate wave function ?7(r) and the field operator 
of uncondcnscd particles ipi{r), satisfying the Bose commutation relations (35). 

The field operator ipl generates the Fock space J-'^ipi), which is orthogonal to the Fock 
space ^(■0) generated by the operator ■^^ The Bogolubov shift (34) reahzes the unitary 
nonequivalent operator representations [32]. Consequently, any operator A[ip], expressed 
through the field operator tp, is defined on the space while the corresponding operator 

74[?7,-?/'i] acts on the space J^^ipi). 

Accomplishing the Bogolubov shift (34), we pass from the Hamiltonian H[iIj] to the 
Hamiltonian H[r],ijji\. Then the average 

<i[7?,V^i] >^ = IVp(r7)i[7y,Vi] (37) 
is defined as the trace over with the statistical operator 

Using average (37), one introduces the quantum-number conservation condition 

<^i(r)>, = 0, (39) 

which defines the operator average 

< 7P{r) >, = ri{r) (40) 

as the order parameter of the system. The condensate wave function is normalized to the 
number of condensed particles 

No= [ \r]{r)\Ur. (41) 



If the system is absolutely stable and the number A^o satisfies the condensation criteria (8) 
or (12), then there exists BEC. 

The Hamiltonian (27), with the infinitesimal term, transforms to H^[r],ilJi\. The related 
statistical operator becomes 

PeW = 7^ r nrr r TIT ' 

Trexp{-/3i/e [77,^/^1]} 
where the trace is over ^(V'l)- Similarly to Eq. (37), we now have the average 

< A[7], ^i] = Tr Pe{v)A[r], ^i] . (43) 

Since 

hm H,[ri,'il:i] = -f/[r;,'0i] , 
averages (37) and (43) are connected through the relation 

lim < A[r), ^i] = < A[r], ^j^] >^ . (44) 
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The Bogolubov shift (34), because of condition (40), exphcitly breaks gauge symmetry. 
When the infinitesimal term is added, then < "01 (r) is not zero. However, in comphance 
with Eq. (44), we have 

lim < i)i{v) = < t/'i(r) >^ = . 

Therefore 

lim < ^(r) >,r, = < ip{r) >n = vi^) , (45) 

which means that gauge symmetry is broken, provided ri{r) is not zero. Thus, the Bogolubov 
shift (34) breaks the gauge symmetry without introducing infinitesimal sources. 

4 Sufficient and necessary conditions 

All observable and thermodynamic quantities of a statistical system can be expressed through 
some correlation functions. Let us consider the class of correlation functions 

C,{il:o,^i) = < . . .4 - ■ -^l ■ -^0- ■ -^l - ■ ■ >e (46) 

consisting of the averages (29) of a product of the field operators -00, ipi, V'o, and ipi in any 
order. And let us introduce the class of correlation functions 

C{ri,^,) = <...rj*...4---V---^i>v (47) 

represented by averages (37). The order of the factors here is the same as in Eq. (46) with 
the change of ipo hy rj. Also, we define the class of correlation functions 

Ceirj^i/Ji) = < ...ri* ...iljl...r]...'ilJi (48) 

given by the averages (43) of a product of rj*, ipl, r], and ipi, whose order is the same as in 
Eq. (47). 

The class of correlation functions (46) is defined for a system, where gauge symmetry 
is broken by means of the infinitesimal source. In this method, the limit £ ^ and the 
thermodynamic limit do not commute, so that 

[lim, lim]C,(V'o,V'i)7^0. (49) 

£— >0 Af— »oo 

The class of correlation functions (48) is constructed by employing the Bogolubov shift (34). 
In the latter method, gauge symmetry is broken explicitly, because of which the limits £ — > 
and N ^ oo now commute, 

[lim lim]C,(77,V'i)=0. (50) 

£— »U Af— >oo 

4.1 Sufficient condition for condensation 

Bogolubov [28,33,34] showed that the methods of breaking gauge symmetry by means of 
infinitesimal sources and by the operator shift (34) are asymptotically equivalent and the 
operator shift (34) is asymptotically exact in the following sense. Let us compare the classes 
defined in Eqs. (46), (47), and (48). The correlation functions from these classes will be 
called similar if they have the same operator structure up to the change of ipo by rj. 
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Bogolubov theorem. 

Similar correlation functions from the classes (46) and (48), in thermodynamic limit, 
coincide, 

lim Ce(V'o,V'i) = lim C,{r|,^^;^) . (51) 

iV— »oo N—*oo 

This limiting equality holds for any e, including £ — > 0. Therefore, taking into account the 
commutation property (50), one has 

lim lim Csiipojipi) = lim C(?7,-?/'i) . (52) 

£— >0 JV— >oo W— >oo 

The limits in the left-hand side of Eq. (52) do not commute. 
From the Bogolubov theorem, it follows, in particular, that 

lim |< Vo(r) >£|' = lim < z/'(l(r)V^o(r) >e ■ (53) 

N—>-oo N^oo 

This is evident, since Eq. (53), with property (51), reduces to the trivial identity \rj\'^ — \r)\'^. 
For a uniform system, Eq. (53) becomes 

lim ^\<ao>e\'^ ^ lim ^ < ajao >e ■ (54) 

This again immediately follows from the Bogolubov theorem (51) and, according to Eq. (33), 
reduces to the identity po = po- Equality (54) has also been obtained by Belyaev [35] and 
Ginibre [36] and recently an elegant proof was given by Lieb et al. [37-39]. 

Equahties (53) and (54) are valid for any £, including e ^ 0. Defining the local condensate 
density 

Po(r) = lim lim < Vo(r)V'o(r) >e , (55) 

£— +U AT— >oo 

we find from Eq. (53) 

lim lim I < Vo(r) >e P = Po(r) . (56) 

The latter equation shows us that if gauge symmetry is locally spontaneously broken, so 
that the left-hand side of Eq. (56) is not identically zero, then there exists local BEG in the 
sense that po(r) > at least for some r. And if there is no spontaneous symmetry breaking, 
at least locally, so that the left-hand side of Eq. (56) is identically zero, then there is no 
BEG, since Po(r) = 0. When spontaneous breaking of gauge symmetry occurs globally, in 
the sense that 

lim lim - / I < ^o(r) >, ^ > , (57) 
then there exists BEG in the usual sense of criteria (8) or (12), since then 



lim — / Pciir) dr > 

V—^0O 

For a uniform system, as is clear, the local and global BEG conditions are equivalent. 

Note that breaking of gauge symmetry can interchangeably be considered for the averages 
of either -0 or ■0o, since ip — ipo + ipi and from the Bogolubov theorem we have 

lim lim < -01(0 >e = lim lim < -01 (r) = lim < -01(0 >r, = , 

£^0 N—i-oo £— >0 Af— >oo N—>-oo 
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where condition (39) is also taken into account. Therefore 



hm hm < ibi^) >e = hm hm < ibnir) >^ . (58) 

The asymptotic exactness of the Bogolubov shift (34) has also been proved for thermody- 
namic potentials by Ginibre [36] and Lieb et al. [37-39] . Let us consider the thermodynamic 
potentials 

= -Tin Trexp{-(3H,[ilj]} , (59) 

with the trace over J-{ip), and 

n,{r]) = -Tin Trexp{-/3i/,[r7,Vi]} , (60) 
with the trace over T{ipi). And let r] be the minimizer for the potential 

Then, for sufficiently general conditions on the interaction potential (see details in Refs. 
[36-39]), one has the Ginibre theorem 

hm 1 [a(77) - fi,] = . (61) 

v — »oo V 

Ginibre [36] proved this theorem for e > as well as for e = 0. 

The Bogolubov and Ginibre theorems show that two ways of breaking gauge symmetry, by 
means of infinitesimal sources and by the Bogolubov shift (34), are asymptotically equivalent 
and that the Bogolubov shift (34) is asymptotically exact in the sense of Eqs. (51), (52), and 
(61). It also follows that the spontaneous gauge symmetry breaking is a sufficient condition 
for the existence of BEG. 

4.2 Necessary condition for condensation 

The fundamental question is whether the gauge symmetry breaking would also be the nec- 
essary condition for BEG. This problem has been discussed in a number of works, of which 
we mention here just a few [1,2,22,26,40]. One often states that gauge symmetry breaking is 
not necessary for BEG and that the letter can occur without any symmetry breaking. This, 
however, is a delusion, though a widespread one. 

Suppose, we do not break gauge symmetry and work with the statistical operator 

^~ Trexp{-PH[7P]} ' ^ > 

in which H[%1)\ is a gauge-symmetric Hamiltonian. The related gauge-symmetric averages 
are defined as 

<A>= Trpi . (63) 
For a uniform Bose-condensed system, Roepstorff [41] proved the inequahty 
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In his proof, he used the assumption that, in thermodynamic hmit, H[^] commutes with -00- 
This is a plausible assumption, since, as is discussed above, in thermodynamic limit, ipQ can 
be replaced by a nonopcrator function rj. A more general proof of inequality (64) was given 
by Lieb et al. in a series of papers [37-39]. Actually, the main requirement, used by Lieb 
et al. [37-39] , was the existence of the Fourier transform for the interaction potential. The 
case of hard-core potentials can also be taken care of by cutting off the potential at some 
finite value that is taken to infinity at the end of the calculations [37-39]. 

In the coordinate representation, the Roepstorff inequality (64) can be written as 

lim - / < Vo(r)V'o(r) >dr < lim lim - / | < ^o(r) >e ?dr . (65) 

I/— >oo V Jv £— »0 V—*oo V Jv 

When the left-hand sides of Eqs. (64) or (65) are nonzero, this means the existence of BEC, 
in compliance with the BEC criteria (8) or (12). But then the right-hand sides of Eqs. 
(64) and (65) are also nonzero, which implies spontaneous symmetry breaking. That is, the 
existence of BEC necessarily assumes spontaneous gauge symmetry breaking. 
Equahties (64) and (65) are the direct analogs of the Griffiths [42] inequahty 

hm < > < hm hm I < S >£ p 

N-*oo s—*0 N—*oo 

for magnetic systems, where 

1 ^ . 

is the mean spin operator and the symmetry breaking term is eNS^, so that e plays the 
role of an external magnetic field. The Griffiths inequahty tells us that the appearance of 
magnetization necessarily implies the spontaneously broken spin-rotational symmetry 0(3). 

In this way, from the Bogolubov theorem, formulated in Eq. (51), equalities (53) to 
(56) follow, which have also been proved by Ginibre [36] and by Lieb et al. [37-39]. These 
equalities show that the spontaneous gauge symmetry breaking is a sufficient condition for 
the occurrence of BEC. And inequalities (64) and (65), proved by Roepstorff [41] and by 
Lieb et al. [37-39], establish that the existence of BEC is necessarily accompanied by the 
spontaneous gauge symmetry breaking. Thence, we come to the conclusion: 

The spontaneous gauge symmetry breaking is the necessary and sufficient condition for 
the existence of Bose-Einstein condensate. 



4.3 Illustration by mean-field model 

In order to illustrate the above conclusion explicitly, let us consider a simple mean-field type 
model with the Hamiltonian 

H = Y^{.^k - l^)aiak , (66) 

k 

where and al are Bose operators and Uk is an effective spectrum, such that (^fe — > 
for all multi-indices k. For a uniform system, the index k becomes the momentum k. But, 
in general, the system can be of any nature. 
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For temperatures T above the condensation temperature Tc, there is no need to break 
gauge symmetry. Then the particle distribution is 

<alak>= [e'^('^'=-'^) (67) 

for any k. The appearance of a Bosc condensate presupposes that the ground-state level ujq 
gets macroscopically populated, as a result of which the number of particles on this level, 



< OoOo > 



-1 



(6^ 



becomes proportional to the total number of particles N . In order that < a\aQ >oc N ^ it is 
necessary that 

lim (cuo - /x) ^ +0 (T < Tc) . (69) 

AT— >oo 

Hence, Eq. (69) is a necessary condition for BEC. 

Breaking gauge symmetry by the infinitesimal term, as in Eq. (27), we now have 

H, = H^e^/N{al + ao) • (70) 

Hamiltonian (70), with H from Eq. (66), can be diagonalized by means of the canonical 
transformation 

ao^bo + — . (71) 

II-CUo 

Then Eq. (70) transforms into 

He^ {ujQ- Mbo + Y^i^k- fJ')alak^ • (72) 

k^Q - ^0 

The latter Hamiltonian is gauge-invariant with respect to the change bo — > 6oe*" and afc — > 
afeC*". Therefore, 

< &o >. = , (73) 
where the averaging is with Hamiltonian (72) and, similarly, 

<ak>e = {k ^ 0) . (74) 

Prom Eq. (71), we find 

< ao >e = . 75) 

For the bilinear averages, we get 

< blbo >e = < ajao > , (76) 
whose form is given by Eq. (68), and 

< alttk >e^< aluk > (ky^O) , (77) 

defined by expression (67). 
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If there is no BEC, hence, condition (69) does not hold, then there is no spontaneous 
gauge symmetry breaking, since 

lini lim ^ =0 il^y^ujo) . (78) 

The number of condensed particles is, generally, different being calculated using the averages 
with gauge symmetry breaking or without it. Let us denote the former as 

No{e) = < ajao >, {e ^ 0) (79) 

and the latter as 

No{0) = < alao > {e = 0) . (80) 
In agreement with transformation (71), we have 

No{s) = < blbo >e +7 ^ . (81) 

From here, keeping in mind equality (76), we get 

< 4ao > = A^o(O) < N^{e) = < ajao >. . (82) 

Let below Tc the condensate arise, so that A''o(£:) > and condition (69) holds. Then, for 
large N, from Eq. (81), we obtain 



T 



2iVo(£ 

For £ 7^ and N ^ oo, Eq. (83) simphfies to 



1 + ^1 + {2(3eYN^{e)N 



(83) 



H^UQ-e \ . (84) 



Noie) ■ 



Using this in transformation (71), we have 



ao ~ 6o - \/No{e) . (85) 
Then, in view of Eq. (73), we obtain 



<ao>e^ -^jNoie) . (86) 
Equation (85) is analogous to the Bogolubov shift, and Eq. (86) yields 

lim < °° ^' I' = l.m ^ . (87) 
The latter equality, together with definition (79), takes the form 

lim ' ^ "° I' = lim 1^2^ , (88) 
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which is equivalent to the Ginibre-Lieb equahty (54). Equation (88) tells us that asymptot- 
ically 

< alao ~ I < ao >e P (A^ oo) . (89) 

Taking in Eq. (88) the limit e — 0, we sec that the spontaneous gauge symmetry breaking 
implies the existence of BEG, being a sufficient condition for the latter. 

On the other hand, from Eqs. (82) and (88), with definitions (79) and (80), we obtain 

hm 1^0^ < lini ' ^ "° I' . (90) 

For e ^ 0, Eq. (90) becomes the Roepstorff-Lieb inequality (64), which establishes the 
necessity of spontaneous gauge symmetry breaking for the BEG existence. 



5 Conservation of the number of particles 

If the spontaneous gauge symmetry breaking is a necessary condition for BEG, then where is 
a flaw in the standard argument of those who claim that BEG docs not need gauge symmetry 
breaking and can be realized without the latter? This argument is as follows. When gauge 
symmetry has been broken, the number-of-particle operator TV does not commute with the 
Hamiltonian H^, hence, N is not an integral of motion. However, in experiments, for instance 
with trapped atoms, the number of particles N can be well controlled and kept constant 
during all the process of measurement. 



5.1 Conservation in the sense of quasiaverages 

The flaw in this argument is that it confuses the operator N with the observable quantity, 
which is the average N =< N >. There is no any contradiction between the noncommu- 
tativity of an operator with a Hamiltonian, in which a symmetry is broken by infinitesimal 
sources, and the fact that the average of this operator can be well controlled and kept con- 
stant in experiments. The situation for BEG is analogous to that for magnetic systems. In 
the latter, the appearance of magnetization is accompanied by the spontaneous breaking of 
spin-rotational symmetry. When the latter is broken, the spin operator S docs not commute 
with the Hamiltonian. This, however, in no way forbids the magnetization, proportional to 
the spin average < S >, from being well controlled, accurately measured, and kept constant 
in experiments. 

To demonstrate what is said above mathematically, let us consider the number-of-particle 
operator 

N = J ij\r)i/j{r) dr . (91) 

This operator commutes with the gauge-symmetric Hamiltonian H[ip], 

[Hm, N] = 0. (92) 

But after breaking the gauge symmetry with the infinitesimal source, the Hamiltonian H^lijj] 
in Eq. (27) does not commute with N. Keeping in mind the separation of the field operator 
into two parts, as in Eq. (17), we have 

[^o(r), n(r')] = V^o(r)v^S(r')^(r') , 
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where n{r) = ■0^(r)'0(r). Prom here, 

[^o(r), ^] = ^o(r). 
Therefore the commutator of = H^lip] with operator (91) is 



[H,, N]=e^pj [^i(r)+V^o(r) 



dY . 



(93) 



The number-of-particle operator (91) is now not an integral of motion, since the right-hand 
side of the Heisenberg equation 

^ I iV = [TV, H,\ (94) 

is not zero, being given by Eq. (93). But the observable quantities are the averages of 
operators. Hence, we should consider the averaged Eq. (94), which gives 



d 
di 



In view of Eq. (93), we get 

d 



dt 



<N >, = ie^ / < ^'^(r) - ^o(r) >e dr . 



(95) 



(96) 



To correctly define thermodynamic limit, we have to work with reduced quantities. For 
example, Eq. (95) can be rewritten as 



lim = t hm 

V^oo dt V V^oo V 

Using Eq. (93) and the Bogolubov theorem (51), we find 



(97) 



d <N>, 
lim 

V^oo dt V 



— ie lim / < rj*{r) — rj{r) >^ dr. 

V-^oo V Jv 



(98) 



Breaking gauge symmetry by infinitesimal sources requires to define all observable quantities 
as Bogolubov quasiaverages. For the condensate function 77 (r), rememebring normalization 
(41) and using the Cauchy-Schwarz inequality, we have 



dr 



< NnV . 



Because of this, 



V Iv'^ ^*(r) - r7(r) dr 
Therefore, from Eq. (98), we obtain 



lim lim — ^ = . 

e^O V^oo dt V 



(99) 
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This means that the average number of particles N, defined by the relation 



N < N >s 
lim — = lim lim , (100) 

does not change in time and, being fixed at the initial time, can be kept constant during any 
experiment. The number of particles N, according to Eq. (100), can be thought as being 
defined by the asymptotic form A'^ ~< N >£, assuming that e ^ and N — > oo. Thus, 
though the operator does not commute with H^, the observable number of particles N is 
a well controlled quantity that can be fixed with any desired accuracy. 

The form of Eq. (98) hints that the influence of the gauge-symmetry breaking source is 
in the appearance of a phase in < TV >£, proportional to e. This is easy to illustrate for a 
uniform system, when 77(r) = ^/po- Then asymptotically, as — > oo, we have 

<[H„ N]>,r^2e^N , (101) 

where 

is the condensate fraction. Keeping in mind that, in accordance with Eq. (100), asymptoti- 
cally N ~< N >£, we obtain from Eq. (98) the asymptotic equation 

dN ^ . X 

— - 2e^ N . 103 
at 

Prom here 

N{t) ~ A^(0) exp {2e^ t) . (104) 
Passing again to quasiaverages, we get 

lim lim = 1 . (105) 

This confirms once more that the observable number of particles can always be treated as a 
well defined, conserved, and precisely controlled quantity. 



5.2 Conservation under Bogolubov shift 

Moreover, the number-of-particle operator can be explicitly an integral of motion, at the 
same time with gauge symmetry being broken. This can be realized in the representation 
employing the Bogolubov shift. The initial Hamiltonian H[il)], which is a functional of the 
field operators ip and ip^ , is gauge invariant, according to property (26). Accomplishing 
the Bogolubov shift (34) leaves the operator structure of the new Hamiltonian H[ijj] = 
H[ri,ipi] the same as that of H[iIj]. Hence, the Hamiltonian H[ri,ipi] is invariant under the 
transformation tjj — > T^e*", that is, rj — > r^e*" and ■01 — V'lC'"- As far as -ff[V'], by assumption 
(26), is gauge invariant, it commutes with the number-of-particle operator iV[?/']. Since the 
Bogolubov shift (34) does not change the operator structure, except for replacing ijj hy ip, 
the Hamiltonian Hlrj, ipi] commutes with the number-of-particle operator 

7V[77,V'i]= U\r)i;{r)dr, (106) 
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so that 

'N[r|,^^], H[ri,ij,]\^0. (107) 
Operator (106) can be written as the sum 

N[ri,^i]=No + Ni (108) 

of the number of condensed particles (41) and the operator 

Ni = J V^I(r)V^i(r) dr (109) 

for the number of uncondensed particles. Thus the number-of-particle operator (106), in view 
of commutator (107), is the integral of motion, which means that the number of particles 

N = <N[rj,iPi]>,, (110) 

with the average defined in Eq. (37), is a conserved quantity, which can be measured and 
fixed with the desired accuracy. At the same time, the gauge symmetry, in compliance with 
Eq. (40), is broken. 

In order to understand how it may happen that the number-of-particle operator is the 
integral of motion, but the gauge symmetry is broken, we should remember that symmetry 
can be broken not solely in the Hamiltonian, by adding infinitesimal sources, but also by 
constructing an operator representation on a space of microstates with broken symmetry 
[31]. To be more precise, we can invoke the decomposition theory [43-45]. 

For a finite system, the space of microstates is the Fock space J-'{4') generated by the field 
operator ■^^ The space is gauge invariant, in the sense that ^^(-06*") = ^("0). The Fock 
space J^{il>i), generated by ipl, is not gauge invariant. Actually, one can define an infinite 
number of spaces 

J^a = :F (V'le^") (0 < a < 27r) . (Ill) 
In thermodynamic limit, the sole space disintegrates into the direct integral 

^ , ^ da 
T 



of the partial spaces (111). 

On each of subspaces (111), one can define a representation of the operator algebra, 
whose members are denoted as 



An = A 



Tie'", V'le*" . (113) 



For this representation, we can construct the corresponding statistical averages by introduc- 
ing the statistical operator 

PocV)) = 7^ 1 nrr ^ ' (H^ 

exp{-l3Ha) 

in which Ha = H [r/e*", ■^le*"], according to Eq. (113). Then we have 

< Aa >r,a = Tr:f Ja{r])Aa . (115) 
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It is easy to notice that 

lim < = < A[r], ipi] >r, , 

a— »0 

with the right-hand side given by Eq. (37). The average of an operator A, defined on the 
gauge-invariant space (112), is decomposed into the integral 

<A>:p= <A^>r,a tt (116) 

Jo ZTT 

over gauge- noninvariant averages defined for subspaces (111). 
In conformity with condition (39), we require that 

<^i(r)>,a = 0. (117) 

Then the gauge-invariant average of -0 is 

/■27r rjfy 

< ^(r) >^ ^ / r;(r)e- — = , (118) 

Jo ZTT 

as it should be. However, the average of the representation ■0a(r) = ■0(r)e*" over the gauge- 
noninvariant space (111) is not zero, 

< Mr) >^a = r?(r)e^" , (119) 

which implies gauge symmetry breaking. This is the way how the gauge symmetry of a 
system can be broken, despite that the number-of-particle operator remains the integral of 
motion. 

One may also notice that, under the assumption that tends, in thermodynamic 

limit, to J^, given by decomposition (112), then 

hm <4{r)Mr)>^\v{r)\' . (120) 

TV— »oo 

This again proves that BEC is equivalent to gauge symmetry breaking. In other words, the 
latter is a necessary condition for the existence of BEC. 



6 Negligible condensate fluctuations 

In recent years much attention has been given to the study of condensate fiuctuations (see re- 
view article [5]). One often claims that these fiuctuations are anomalous or even catastrophic. 
However, if to calculate these fluctuations correctly, not forgetting about the broken gauge 
symmetry, which is a necessary condition, then they are neither anomalous nor catastrophic, 
but, moreover, they are negligible in thermodynamic limit. 

One describes the condensate fiuctuations as the fiuctuations associated with the con- 
densate number-of-particle operator 

No = J 4ir)M^) dr . (121) 

The fiuctuations are characterized by the dispersion of operator (121), 

Al{No) = <N^>e-< No >l . (122) 
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The dispersion of an extensive operator is not, of course, an observable quantity by its own. 
The observable quantity is the ratio 

Un. ^ . lim lim ^ . (123) 

N-*oo N £-*0 N->-oo N 

which is analogous to the definition of other extensive observables [46,47]. Denoting 

no(r) = 4{r)Mr) , (124) 
we can write down dispersion (122) as 

AliNo) = 1 [< no(r)no(r') >, - < no(r) >,< no(r') >,] drdr' . (125) 

In thermodynamic limit, according to the Bogolubov theorem, we should make the replace- 
ment iPo{t^) ^(i")) hence 

no{r) po(r) = |77(r)|^ . 

Then Eq. (125) immediately gives A^(iVo) ~ 0. 

One may ask what happens if one represents dispersion (125) is a slightly different form 
following from the commutation relation (19)? Using the latter, we can write 

no(r)no(r') = ^i(r)4(r')^o(r')^o(r) + 4ir)M^')M^)iplir') . (126) 

Wc should expect that the result for dispersion (125) could not depend on the used commu- 
tation relation (19), since the condensate operators ipo{r) and ipQ^r'), in agreement with Eq. 
(23), commute in thermodynamic limit. More accurately, there is the property 



lim < 



^o(r),V5(r') >e=0. (127) 



One should exercise extreme caution integrating expressions containing the commutation 
relations (19). As is explained in Section 2, such expressions have to be understood in 
the sense of equations for distributions. Their integration is correctly defined only for the 
integration with an integrable function, as is specified in Section 2. Thus the integration of 
Eq. (127) over the spatial variables has sense only for the case 

lim / < Uo{r),4{r')] >e /(r, r') dr dr' = , (128) 

AT— >oo Jv L 

where /(r,r') is an integrable function, such as defined in subsection 2.2. The integration 
without such a function can yield wrong results. For example, we could get 



/ <[Mr),4{r') 

J V 



>F dr 



for large N oo, which would contradict the commutativity property (23). To exclude the 
appearance of such incorrect terms, one has to employ the regularized limit 

hm ^ = 

N^oo N 
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lim lim ^ / [< no(r)no(r') >e - < rio(r) >^< no(r') /(r, r') dr c/r' 



(129) 



Reducing the regularizing function /(r,r') to unity can be done only a/ter thermodynamic 
limit, but not before, since 

lim, lim ^ . 

Suppose, we decide to use form (126) in dispersion (125). Employing the equality ■0o(r) = 
ao(/?o(i"), the Bogolubov replacement Oq — > condition (22), and taking care of correct 

definitions, we have 



lim — 



V 



< lim 



<^^(r)^o(rO >. ipo{r)LpQ{r')f{r,r') dr dr' 
const 



< 



Po 



- /(r,r') dvdv' 



. 



Thence, the additional term in Eq. (126), caused by the commutation relations (19), does 
not contribute to the reduced dispersion (129). In any case, we have 



N-^oo N 



(130) 



This means that the condensate fluctuations are negligible in thermodynamic limit. 

If we would employ representation (15) in the definition of operator (121), we would get 



7V„ = /n„(r)<ir = aJ«„. 



Hence, 



Nl = ajooajoo . 

Using here the commutation relation [oq, Oq] = 1, we have another form 



(131) 
(132) 

(133) 



It looks hke Eqs. (132) and (133) give different answers for dispersion (125). For instance, 
keeping in mind the Bogolubov replacement VNq and using form (132), we have 

A^(iVo) = 0, in agreement with Eq. (130). At the same time, using form (133), we would 
get A2(iVo) No, which contradicts Eq. (130). The cause of this contradiction is that the 
spatial integrations, in the process of passing from the spatial representation to Eq. (131), 
have been done without a regularization function. To correct this defect, one has either to 
return to the spatial representation and to proceed as has been explained above or to invoke 
the following argumentation. 

The Bogolubov replacement Oq ^/No implies that the operators ao and should 
commute in thermodynamic limit. More correctly. 



lim ^Kl^^o, 

AT-^oo No 



(134) 
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where A^o =< -^o >£ tends to infinity proportionally to iV, as iV — > oo. To correctly define 
< TVq >£, one has to consider the reduced quantity < Nq >^ /Nq. Then, form (132) gives 



lim 



>e 



^ 1 . 



Form (133) yields the same answer 



lim 

AT— >oo 




N-^oo 



lim 1 + 




In that sense, one always has 



< ^0 




(iV ^ oo) , 



(135) 



independently of the form involved. Thus one always returns to Eq. (130) telling us that 
the condensate fluctuations are thermodynamically negligible. 

One could come to catastrophic condensate fluctuations only by means of incorrect cal- 
culations, assuming the existence of BEG but without breaking gauge symmetry. In the 
latter case, using the Wick theorem, one would get 



This what one calls the catastrophic condensate fluctuations, since the normal dispersion 
must be proportional to N. One often blames the grand ensemble for this unreasonable 
behavior of the dispersion. However, as has been just explained, it is not the grand en- 
semble to blame, but this fictitious catastrophe is merely provoked by wrong calculations. 
The correct way of treating the problem results in thermodynamically negligible condensate 
fluctuations in the sense of Eq. (130). It is worth stressing that result (130) is general and 
is valid for arbitrary Bose systems, whether uniform or nonuniform and even whether they 
are equilibrium or nonequilibrium. 

7 Representative statistical ensembles 

The basic point, analyzed and explained in this paper, is that spontaneous gauge symmetry 
breaking is the necessary and sufficient condition for the existence of BEC. Although it 
is, probably, possible to calculate some particular quantities of a Bose-condensed system 
without explicit symmetry breaking, but, in general, such a way of calculation is not correct 
and can lead to wrong results. There are two methods of gauge symmetry breaking, by 
introducing into the Hamiltonian infinitesimal sources or by means of the Bogolubov operator 
shift. These methods are equivalent in thermodynamic limit. However, before the latter, the 



< ATq > = < ajao > (2-F < ajao >) , 



from where 



with A^o =< OoOo >oc A^". As a result, one finds 



< TVq^ > - < iVo >^ oc . 



23 



methods are not equivalent, since they involve different operator variables possessing different 
commutation relations. To formulate this difference in a precise form, it is convenient to 
invoke the notion of representative statistical ensembles [31,46,48-50]. 

A statistical ensemble is a pair {J^, p} of the space of microstates and a statistical 
operator p. The latter, for an equilibrium system, is given by the Gibbs form characterized 
by the appropriate Hamiltonian, and for a nonequilibrium system, p is described by the 
Liouville equation, whose dynamics is governed by the related Hamiltonian. 

A representative statistical ensemble is an ensemble equipped with all additional condi- 
tions that are necessary for a unique description of the given physical system [31,46]. Addi- 
tional conditions may include symmetry properties, normalization conditions, conservation 
conditions, and so on. Actually, only representative ensembles have sense. 

Suppose, we are considering a Bose system in the frame of the field-operator technique. 
With a field operator ip, the space of microstates is the Fock space J-'{ip) generated by ipK 
The field operators, as is well known, do not conserve the number of particles. Therefore, 
working with the field-operator technique, it is necessary to use the grand canonical ensemble 
fixing the average number of particles N =< N >. The statistical operator p is defined by 
a Hamiltonian H[iIj] which is a functional of ip. To emphasize this, we shall write p[H[ijj)]. 
Let us denote by H[ip] the energy operator. And let us define the statistical ensemble 

{J^{^), p{H[m , H[^] = Hm - f^N . (136) 

This is the standard grand canonical ensemble with the grand Hamiltonian H[ijj], which is 
assumed to be gauge invariant. Ensemble (136) is representative for a system without BEG. 
But as soon as the latter appears, ensemble (136) is not anymore representative, since BEG 
necessarily requires gauge symmetry breaking. Accomplishing this by means of the method 
of infinitesimal sources, adding a gauge-symmetry breaking term F, we have the statistical 
ensemble 

m^), pmm , = hhj] -pN+et. (m) 

This ensemble is representative for a system with BEG, the operator averages being defined 
as the Bogolubov quasiaverages. The number of particles is now given by N —< N >^. 
The chemical potential p is the Lagrange multiplier preserving the normalization condition 
for the number of particles A^. In ensemble (137), there is the sole field variable ijj and the 
only normalization condition for N. Hence, just one Lagrange multiplier is required. The 
representation (17) for ip, as a sum of two terms, is rather formal, since neither of these 
terms comprises an independent variable. The parts of sum (17) are interconnected through 
their commutation relations (19), (20), and (21). It is, of course, possible to work with 
ensemble (137), invoking perturbation theory in powers of particle interactions, as has been 
done by Belyaev [35]. This, however, has two weak points. First, such a perturbative way 
makes it possible to deal exclusively with weakly interacting systems. Second, this way is 
very cumbersome because of the complicated commutation relations (19) and (20). 

It is tempting to resort to the Bogolubov theorem and to make the Bogolubov shift 
(34), passing to the simple commutation relation (35). However, one has to remember that 
the Bogolubov shift becomes exact only in thermodynamic limit, but not before it. At the 
same time, all calculations are usually accomplished before thermodynamic limit. Using 
the commutation relations (35) before this limit is incorrect, makes the theory not self- 
consistent, and leads to different paradoxes. The most dangerous such a paradox is the 
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Hohenberg-Martin dilemma of conserving versus gapless theories [51]. This paradox has 
recently been resolved in Rcfs. [32,46,48-50,52,53]. 

The solution to the problem is merely the necessity of being mathematically accurate. 
Yes, we can employ the Bogolubov shift (34) and to use the simple commutation relations 
(35) even before thermodynamic limit, provided N is large. But then it is necessary to 
understand that the Bogolubov shift introduces two independent variables, the condensate 
wave function ^^(r) and the operator of uncondensed particles ipi{r). So, there are now two 
normalization conditions, for the number of condensed particles 

No^J \vir)\Ur (138) 

and for the number of uncondensed particles 

N, = < TVi >^ , (139) 

where Ni is given by Eq. (109) and the average is defined in Eq. (37). Preserving two 
normalization conditions requires two Lagrange multipliers. In addition, one must be aware 
of the fact that the Bogolubov shift (34) induces a new operator representation, which is 
nonequivalent to that used in ensemble (137). Now, the space of microstates is the Fock 
space J^{ipi) generated by V'l- The space J^{ipi) is asymptotically orthogonal to the space 
J^{ip) of ensemble (137), with all operators forming the new representation on this space 
[32]. After accomplishing the Bogolubov shift, we come to the representative ensemble 

p{H[r^,^l^,])} , H[ri,^,] = H[rj,^^] - fXoNo - f^iNi , (140) 

with the statistical operator (38) and the grand Hamiltonian H[r],ilJi\ which is a functional 
of two variables and which contains two Lagrange multipliers, /xq and /ii, guaranteeing 
the validity of two normahzation conditions (138) and (139). Since the Bogolubov shift 
explicitly breaks gauge symmetry, we do not need to invoke infinitesimal sources. The 
representative ensemble (140) is the basis for developing a completely self-consistent theory 
of Bose-condensed systems [46,48-50], which is free of paradoxes and which enables the 
description of strongly interacting systems with BEC [52,53]. 

The case of Bose-condensed systems serves as a very illustrative example emphasizing the 
principal importance of defining a representative statistical ensemble for the system under 
consideration. It may happen that several ensembles can correctly describe a statistical 
system. Then one says that these ensembles are equivalent. For instance, the majority of 
statistical systems can be equivalently treated by means of either grand canonical, canonical, 
or microcanonical ensembles. But this is not compulsory. Even more, pronouncing just 
the term of an ensemble, say, grand canonical, is not sufficient, since, as has been shown 
above, there can exist several grand canonical ensembles. Thus ensembles (136), (137), and 
(140) all are grand canonical, though they are different and not all equivalent. The most 
important is that the employed ensemble be representative. In the other case, it is unable 
to correctly represent the considered system. Several ensembles are equivalent if and only if 
they are representative [46]. There are some models with long-range interactions for which 
the microcanonical and canonical ensembles are not equivalent [54,55]. But this in no way 
is a fact that would shake the basis of statistical mechanics. This just means that one or 
both of the used ensembles are not representative. Say, the microcanonical ensemble may 
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occur to be not representative for such mean-field models. Or it may be that some additional 
constraints are to be imposed in order to correctly define the model and to make ensembles 
representative, hence, equivalent. 

In conclusion, it is worth emphasizing that the principal problems analyzed in the present 
paper concern any kind of Bose systems with BEG, whether uniform or nonuniform. These 
could be uniform Bose gases and superfiuid liquids or nonuniform clouds of trapped Bose 
atoms. As is mentioned in the Introduction, multiquark clusters can also form BEG in the 
interiors of neutrons stars and in fireballs created in heavy-ion collisions [17-19]. Another 
example is the BEG of excitons [56,57]. Nonuniform condensates can appear in a uniform 
system, when the spectrum of collective excitations [58,59] or the single-particle spectrum 
[60-62] touch zero at the finite momentum ko. Then a condensate with a periodic structure, 
prescribed by the momentum kg, arises, forming a coherent crystal [58,59]. A similar periodic 
density evolves in nuclear matter under pion condensation [63-72] . Periodic condensates can 
occur, in general, for any Bose particles in external fields [71,73]. Finiteness of a sample 
also can impose boundary conditions that would induce a periodic BEG, such as happens 
in magnetic films under magnon condensation [74,75]. The abundance of different types of 
BEG makes crucial the necessity of understanding the basic features of its correct theoretical 
description. 
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between BEG and gauge symmetry breaking. 
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